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We investigate the influence of external laser excitations on the average current through bridged 
molecules. For the computation of the current, we use a numerically very efficient formalism that 
is based on the Floquet solutions of the time-dependent molecule Hamiltonian. It is found that the 
current as a function of the laser frequency exhibits characteristic peaks originating from resonant 
excitations of electrons to bridge levels which are unoccupied in the absence of the radiation. The 
electrical current through the molecule can exhibit a drastic enhancement by several orders of 
^1^. magnitude. 
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I. INTRODUCTION 



In a seminal workJ3 Aviram and Ratner proposed almost thirty years ago to build elements of electronic circuits 
— in their case a rectifier — with single molecules. In the present days theixp-vision starts to become reality and 
the experimental and theoretical study of such systems enjoys a vivid activity .Ba'cl Recent experimental progress has 
O ' enabled reproducible measurementaj'cl of weak tunneling currents through molecules which are coupled by chemisorbed 
' thiol groups to the gold surface of external leads. A necessary ingredient for future technological applications will be 
^ the possibility to control the tunneling current through the molecule. 

^ Typical energy scales in molecules are in the optical and the infrared regime, where basically all of the today's 

. lasers operate. Hence, lasers represent a natural possibility to control atoms or molecules and also currents through 

them. It is for example possible to induce by the laser field an oscillating current in the molecule which under certain 
^H ^l asymmetry con<^tions is rectified by the molecule resulting in a directed electron transport even in the absence of any 
O ^ applied voltage.B-B Another theoretically predicted effect is the current suppression by the laser ffelds which offers the 
I— possibility to control and switch the electron transport by light. Since the considered frequencies lie below typical 
plasma frequencies of metals, the laser light will be reflected at the metal surface, i.e. it does not penetrate the leads. 
^ ' Consequently, we do not expect major changes of the leads' bulk properties — in particular each lead remains close 
, to equilibrium. Thus, to a good approximation, it is sufficient to consider the influence of the driving solely in the 
\^ ' molecule Hamiltonian. In addition, the energy of infrared light quanta is by far smaller than the work function of a 
, common metal, which is of the order of 5 eV. This prevents the generation of a photo current, which otherwise would 
Q\ ' dominate the effects discussed below. 

' Recent rt^jeoretical descriptions of the molecular conductivity in non-driven situations are based on a scattering 
i approach, t^ltJ or assume that the underlying transport mechanism is an electron transfer reaction from the donor 
' to the acceptor site and that the conductivity can be derived from the corresponding reaction rates .Q It has been 
^ , demonstrated that both approaches yield essentially identical results in a large parameter regime. L2 -Within a high- 
' temperature Hmit, the electron transport on the wire can be described by inelastic hopping events £3o ■— 1|— ■ 
' Atoms and molecules in strong oscillating fields have been widely studied within a Floquet formalism.E§li3 This 
^ ^H ^' suggests utilizing the tools that have been acquired in that area, thus, developing a transport formalism that combines 
^_ Floquet theorjjJbE^ a driven molecule with the many-particle description of transport through a systerrL±|ia± is coupled 
• • . to ideal leads E'EtlI Such an approach is devised much in the spirit of the Floquet-Markov theorytZlEJ for driven 
. ' dissipative quantum systems. 

I II. FLOQUET TREATMENT OF THE ELECTRON TRANSPORT 

The entire system of the wire in the laser-field, the leads, and the molecule-lead coupling as sketched in Figure |^ is 
described by the Hamiltonian 

-leads ■ (1) 



The wire is modeled by N atomic orbitals |n), n — 1, . . . ,iV, which are in a tight-binding description coupled by 
hopping matrix elements. Then, the corresponding Hamiltonian for the electrons on the wire reads in a second 
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FIG. 1: Level structure of a molecular wire with N = 8 atomic sites which are attached to two leads. 



quantized form 

i^wirc (t) = ' (2) 

n,n' 

where the fermionic operators c„, cj^ annihilate, respectively create, an electron in the atomic orbital \n) and obey the 
anti-commutation relation [c„, cl^i]+ = Sn,n'- The influence of the laser field is given by a periodic time-dependence of 
the on-site energies yielding a single particle Hamiltonian of the structure Hnn'{t) = Hnn'{t + T), where T = 2ti/VI 
is determined by the frequency U, of the laser field. 

The orbitals |1) and \N) at the left and the right end of the molecule, that we shall term donor and acceptor, 
respectively, are coupled to ideal leads (cf. Fig. |^) by the tunneling Hamiltonians 

-ffwiro-loads = '^{'^qL c\i^Ci + VqR c\^Cn) + H.C. (3) 

q 

The operator CqL (cqR) annihilates an electron in state Lq (Rq) on the left (right) lead. The leads are modeled as 
non-interacting electrons with the Hamiltonian 

q 

where e^^ is the single particle energy of the state qL and correspondingly for the right lead. As discussed above, the 
leads can be described by a grand-canonical ensemble of electrons, i.e. by a density matrix 

fillcads.cq « exp [-(i/lcads - ML^L - Mfl^i?,)/^^^] , (5) 

where ^Il/r ^^re the electro-chemical potentials and Ni^/r = ^qC^qi^^j^Cqi^/R the electron numbers in the left/right 
lead. As a consequence, the only non-trivial expectation values of lead operators read 

(clL^^qL) ^ f{<^qL - ^^l)- (6) 

Here, f{x) = (1 + e^/*^^-'")^^ denotes the Fermi function. 



A. Perturbation theory 



While the leads and the wire, including the driving, will be treated exactly, we take the wire-lead Hamiltonian as a 
perturbation into account. Starting from the Liouville-von Neumann equation ifig{t) — [H{t), g{t)] together with the 
factorizing initial condition p(to) = f?wirc(io) ® f?ioads,cq, we derive by standard techniques an approximate equation of 
motion for the total density operator Q{t). This is most conveniently achieved in the interaction picture with respect 
to the uncoupled dynamics where the Liouville-von Neumann equation reads 



(7) 
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The tilde denotes the corresponding interaction picture operators, X(t,t') = U}j{t,t') X{t) Uo{t,t') where the propa- 
gator of the wire and the lead in the absence of the lead-wire coupling is given by the time-ordered product 

Uo{t,t') = Texp f dt" [i?wirc(t") + i^lcads]) • (8) 



Equation is equivalent to the following integral equation 



i /■* 

g{t,tQ) = g{to,to) - - dt'[H^ii.c-icads{t',to),g{t',to)]. (9) 

'to 



n 



We reinsert this expression into the differential equation (0) and use that to zeroth order in the molecule-lead coupling 
the interaction-picture density operator does not change with time, Q{t — t.Iq) w g{t,to)- A transformation back to 
the Schrodinger picture results in the following approximate equation of motion for the total density operatorQa 



i 1 f ~ 

Q(t) = - -[H^irc{t) + -fflcads, g{t)] - /dT[7Jwire-leads, [-ff wire-leads " T,t), g{t)]] 



(10) 



[-^^wirc — Icadsi Uo{t,to)g{to)U^o(t^to)]- 



Since we only consider asymptotic times t — to — > oo, we have set the upper limit in the integral to infinity. The 
third term in Eq. ( p^ ) stems from the initial condition at to in the integrated form (||) of the Liouville-von Neumann 
equation. For the chosen factorizing initial condition, it will not contribute to the expectation values calculated below. 

The net (incoming minus outgoing) current through the left contact is given by the negative time derivative of the 
electron number in the left lead, multiplied by the electron charge — e, i.e. 

lL{t) ^ etT[g{t)NL]. (11) 

We insert g{t) from Eq. ( |l0| ) and obtain an expression that depends on the density of states in the leads times their 
coupling strength to the connected sites. At this stage it is convenient to introduce the spectral density of the lead-wire 
coupling 

^L/Rie) - ^ E \yiL/B.\'Sie - e,i/«), (12) 

which fully describes the leads' influence. If the lead states are dense, ri//j(e) becomes a continuous function. Because 
we are mainly interested in the behavior of the molecule and not in the details of the lead- wire coupling, we assume 
that the conduction band width of the leads is much larger than all remaining relevant energy scales. Consequently, 
we approximate in the so-called wide-band limit the functions Ti:^/ii{e) by the constant values r^/zj. After some 
algebra, we flnd for the time- dependent net electrical current through the left contact the expression 

C30 

lL{t) ^^ReJdTjde e''^/^[{d\{t - r, t) ci) - [ci,d\{t - r, t)]+f{e - ml)}, (13) 



and correspondingly for the current through the contact on the right-hand side. Here, we made the assumption, that 
the leads are at all times well described by the density operator (||). Note that the anti-commutator [ci, c\{t — r, t)]^ 
is in fact a c-number. Like the expectation value {c\{t — r, t) ci) it depends on the dynamics of the isolated wire and 
is influenced by the external driving. The first contribution of the e-integral in Eq. ( p^ ) is readily evaluated to yield 
an expression proportional to 5{t). Thus, this term becomes local in time and reads eri(c|ci). 



B. Floquet decomposition 

Let us next focus on the single-particle dynamics of the driven molecule decoupled from the leads. Since its 
Hamiltonian is periodic in time, Hnn>{t) = Hnn'{t + T), we can solve the corresponding time-dependent Schrodinger 
equation within a Floffafi t . anpij oach. This means that we make use of the fact that there exists a complete set of 
solutions of the formN'taMBy 



|*a(<)) =e~-°*/''|$,(i)), |$„(i)) = Wo.{t + r)) 



(14) 
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with the quasi-energies £„. Since the so-called Floquet modes \^a{t)) obey the time-periodicity of the driving field, 
they can be decomposed into the Fourier series 

=^e-*'="*|<i>„,fe). (15) 

A; 

This implies that the quasienergies Eq come in classes, 

ea,k^ea + khn, 0,±1,±2,..., (16) 

of which all members represent the same solution of the Schrodinger equation. Therefore, the quasienergy spectrum 
can be reduced to a single "Brillouin zone" —hfl/2 < e < fiil/2. In turn, all physical quantities that are computed 
within a Floquet formalism are independent of the choice of a specific class member. Thus, a consistent description 
must obey the so-called class invariance, i.e. it must be invariant under the substitution of one or several Floquet 
states by equivalent ones, 

Co, |$a(i)> ea + ko^nn, e^'=°"'|$„(t)), (17) 



where ki, . . . , kpf are integers. In the Fourier decomposition (|15|), the prefactor exp{ikaflt) corresponds to a shift of 
the side band index so that the class invariance can be expressed equivalently as 



1 I ^ a, 



kahfl, |$a,fe+fcj. (18) 

Floquet states and quasienergies can be obtained from the quasienergy equationlll0lll'00@ 

d 

\n)nnn'(t){n'\ - in 



A wealth of methods for the solution of this eigenvalue problem can be found in the literature. For an overview, we 
refer the reader to the reviews in Refs. ^ 22, and the references therein. 

As the equivalent of the one-particle Floquet states \^a{t)), we define a Floquet picture for the fcrmionic creation 
and annihilation operators cj^, c„, by the time-dependent transformation 

n 

The inverse transformation 

c„ = ^(n|$„(t))c„(t) (21) 



follows from the mutual orthogonality and the completeness of the Floquet states at equal times.ll3c3 Note that the 
right-hand side of Eq. ( pT| ) becomes independent after the summation. The operators Cq (t) are constructed in such 
a way that the time-dependences of the interaction picture operators Ca (t — r, t) separate, which will turn out to be 
crucial for the further analysis. Indeed, one can easily verify the relation 

icit -T,t)^ U^{t - T, t) cc.it - t) Uo{t - r, t) 
= e"°^/''c„(t) 

by differentiating the definition in the first line with respect to r and using that |$Q.(i)) is a solution of the eigenvalue 
equation (|9|). The fact that the initial condition Ca{t,t) = Ca{t) is fulfilled completes the proof. The correspon ding 
expression for the interaction picture operator in the on-site basis, c„(t — T,t), can be derived with help of Eq. ( |2l| ) 
at time t — t together with ( |2^ ) to read 

5„(i -T,t)^ Y.^n\-f^it - r))e-°^/''c„(i) (23) 

a 

= ^e*(^"/'"'+'=")^e-*'^'"*(n|$„,fc)c„(t). (24) 

ak 

Equation s (p^ ), (p^), consequently allow to express the interaction picture operator c\(t — T, t) appearing in the current 
formula (|13[) via Ca{t), dressed by exponential prefactors. 
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This spectral decomposition allows one to carry out the time and energy integrals in the expression 
current entering the wire from the left lead. Thus, we obtain 



L3h for the net 



(25) 



with the corresponding Fourier components 



If = er. 



i^ck' + k" 1 1) (l|«'/3,fc+fc")^a/3,fc' 

afSk'k" 

) + (^a,A.-'-A;|l)(l|^Q,A;'))/(^Q,fc' ^ f^L) 



Oik' 



(26) 



Here, we have introduced the expectation values 

Ro.p{t) ^ {ci{t)cp{t)) = Rl^{t) 

E-iknt p 
e riap,k- 



(27) 
(28) 



The Fourier decomposition in the last line is possible because all Rapit) are expectation values of a linear, dissipative, 
periodically driven system and therefore share in the long-time limit the time-periodicity of the driving field. In the 
subspace of a single electron, i?^^ reduces to the density matrix in the basis of the Floquet states which has been 
used to describe dissipative driven quantum systems in Refs. |l^,|l^,|3 24 25|,p6|. 

The next step towards the stationary current is to find the Fourier coefficients Rap.k at asymptotic times. To this 
end, we derive from the equation of motion ( [To| ) a master equation for Rap (t) . Since all coefficients of this master 
equation, as well as its asymptotic solution, are T-periodic, we can split it into its Fourier components. Finally, we 
obtain for the Rap.k the inhomogeneous set of equations 



■a(3',k' 



(29) 



IB'k" 



+ X] (*a',fc'+fc" |1> {M^a,k+k")Ra'f3.k' 
a'k" 

- ($/3,fc'-fe|l)(l|$a,fc'>/(ea,fe' - Mi) 

- ($^^fc'|l)(l|$a.fc'+fc)/(e/3,fc' - Mi) 

-|- same terms with the replacement 

{ri,^i,|i>(i|} ^ {^R,^^R,\N){N\}. 

For a consistent Floquet description, the current formula together with the master equation must obey class invariance. 
Indeed, the simultaneous transformation with (18) of both the master equation (p9[) and the current formula (|2l) 



amounts to a mere shift of summation indices and, thus, leaves the current as a physical quantity unchanged. 

For the typical parameter values used below, a large number of sidebands contributes significantly to the Fourier 
decomposition of the Floquet modes \^a{t)). Numerical convergence for the solution of the master equation ([2S|), 
however, is already obtained by just using a few sidebands for the decomposition of Raoit). This keeps the numerical 
effort relatively small and justifies a posteriori the use of the Floquet representation (pl|). Yet we are able to treat 
the problem beyond the rotating-wave-approximation. 



C. Time-averaged current through the molecular wire 

Equation ( p5| ) implies that the current /l(0 obeys the time-periodicity of the driving field. Since we consider here 
excitations by a laser field, the corresponding driving frequency lies in the optical or infrared spectral range. In an 
experiment one will thus only be able to measure the time-average of the current. For the net current entering through 
the left contact it is given by 

Il = II = cTlY, [Y,{<^a,k'+k\l)m0,k')Ro.p,k - (1>a,fe|l)(l|<I>a,fe)/(ec.,fe - A^l)] . (30) 

ak pk' 
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By replacing {|1), _L} {\N), R}, one obtains for the current which enters from the right, iRit), and the corresponding 
Fourier coefficients and time averages. 

Total charge conservation of the original wire-lead Hamiltonian (|l|) of course requires that the charge on the wire 
can only change by current flow, amounting to the continuity equation (5wire(0 = I hit) + lR{t). Since asymptotically, 
the charge on the wire obeys at most the periodic time-dependence of the driving field, the time-average of Owire (t) 
must vanish in the long-time limit. From the continuity equation one then finds that II + Ir = 0, and we can 
introduce the time-averaged current 

I^Il^-Ir. (31) 



This continuity equation can be obtained directly from the average current formula (30) together with the master 
equation (p9|), as has been explicitly shown in Ref. 0. 



III. LASER-ENHANCED CURRENT 



A. Bridged molecular wire 

As a working model we consider a molecule consisting of a donor and an acceptor site and N — 2 sites in between 
(cf. Fig. |l|). Each of the sites is coupled to its nearest neighbors by a hopping matrix elements A. The laser field 
renders each level oscillating in time with a position-dependent amplitude. Thus, the corresponding time-dependent 
wire Hamiltonian reads 

Hnn'{t) = - A{6n,n' + 1 + Sn+l,n') + [En ^ Ax^ COs(m)] Sn„' , (32) 

where a;„ — {N + 1 — 2n)/2 is the scaled position of site \n). The energy a{t) equals the electron charge multiplied 
by the electrical field amplitude of the laser and the distance between two neighboring sites. The energies of the 
donor and the acceptor orbitals are assumed to be at the level of the chemical potentials of the attached leads, 
El = En = hl — Hr- The bridge levels E„ , n — 2, . . . , N — 1, lie Eb ^ A above the chemical potential, as sketched 
in Figure |l|. 

In all numerical studies, we will use a symmetric coupling, F^ = Tr = F. The hopping matrix element A serves 
as the energy unit; in a realistic wire molecule, A is of the order 0.1 eV. Thus, our chosen wire-lead hopping rate 
F = 0.1A/?i yields eF = 2.56 x 10~^ Ampere and Vl « 10A/?i corresponds to a laser frequency in the near infrared. 
For a typical distance of SA between two neighboring sites, a driving amplitude A = A is equivalent to an electrical 
field strength of 2 x 10^ V/cm. 



B. Average current at resonant excitations 



Let us first discuss the static problem in the absence of the field, i.e. for A = Q. In the present case where the 
coupling between two neighboring sites is much weaker than the bridge energy, A <C one finds two types of 
eigenstates: One forms a doublet whose states are approximately given by (|1) ± \N))/^/2. Its splitting can be 
estimated in a perturbational approachEO and is approximately given by 2A(A/i?B)^^^. A second group of states is 
located on the bridge. It consists of iV — 2 levels with energies in the range [Er — 2A, Er + 2A]. In the absence of 
the driving field, these bridge states mediate the super-exchange-between the donor and the acceptor. This yields an 
exponentially decaying length dependence of the conductanceElld 

This behavior changes significantly when a driving field with a frequency k, Er /h is switched on. Then the 
resonant bridge levels merge with the donor and the acceptor state to form a Floquet state. This opens a direct 
channel for the transport resulting in an enhancement of the electron current as depicted in Figure ^ where we plot 
the current amplification, defined as the ratio of the time-averaged current to the current in the absence of the laser, 
•q = I / Iq- In a wire with N = 8 sites, one finds peaks in the current when the driving frequency matches the energy 
difference between the donor/acceptor doublet and one of the — 2 = 6 bridge levels. The applied voltage is always 
chosen so small that the bridge levels lie below the chemical potentials of the leads. The amplification, can assume 
many orders of magnitude, cf. Figure |^. Generally, the response of a system to a weak resonant driving scales with the 
damping and the driving amplitude. Figure]^ demonstrates this behavior for the peaks of the electrical current. The 
peak heights at the maxima of the time-averaged current are found proportional to A^/F. A further scaling behavior is 
found for the current peaks as a function of the wire length: The average current no longer exhibits the exponentially 
decaying length dependence that has been found for bridged super-exchange. By contrast, it emerges proportional to 
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FIG. 2: Amplification of the time-averaged current through the wire sketched in Figure |l]with Eb = lOA. The scaled amplitude 
is A = O.IA; the applied voltage fiL — ftR = 5A/e. The other parameters read F = Fl = F_r — O.lA/h, ksT — 0.25A. 
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FIG. 3: Average current / as a function of the the driving frequency fl for various driving amplitudes A and coupling 
strength F = Fz, = Fij. All the other parameters are as in Fig. 0. 
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FIG. 4: Average current / as a function of the the driving frequency SI for various wire length A''. All the other parameters 
are as in Fig. 0. 
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1/(7V — 1). This can be appreciated in Figure Q where the scale of the abscissa is chosen proportional to — 1 such 
that it suggests a common envelope function. Put differently, the current is essentially inversely proportional to the 
length as in the case of Ohmic conductance. 

In summary, we find current peaks whose height /peak scales according to 

^2 

/peak OC (33) 

Thus, the current is especially for long wires much larger than the corresponding current in the absence of the driving. 



IV. CONCLUSIONS 

We have presented a detailed derivation of the Floquct transport formalism which has been applied in Refs. 
The analysis of a bridged molecular wire revealed that resonant excitations from the levels that connect the molecule 
to the external leads to bridge levels yield peaks in the current as a function of the driving frequency. In a regime 
with weak driving and weak electron-lead coupling, A 3> F, ^, the peak heights scale with the coupling strength, 
the driving amplitude, and the wire length. The laser irradiation induces a large current enhancement of several 
orders of magnitude. The observation of these resonanceSji»uld serve as an experimental starting point for the more 
challenging attempt of measuring quantum ratchet effectsEra or current switching by laser fields.a 



V. ACKNOWLEDGEMENT 



This work has been supported by SFB 486 and by the Volkswagen-Stiftung under grant No. 1/77 217. One of us 
(S.C.) has been supported by a European Community Marie Curie Fellowship. 



^ A. Aviram and M. A. Ratner, Chem. Phys. Lett. 29, 277 (1974). 

^ C. Joachim, J. K. Gimzewski, and A. Aviram, Nature 408, 541 (2000). 

3 A. Nitzan, Annu. Rev. Phys. Chem. 52, 681 (2001). 

* P. Hanggi, M. Ratner, and S. Yaliraki, Special Issue: Processes in Molecular Wires, Chem. Phys. 281, pp. 111-502 (2002). 
^ X. D. Cui et al, Science 294, 571 (2001). 

^ J. Reichert, R. Ochs, D. Beckmann, H. B. Weber, M. Mayor, and H. v. Lohneysen, Phys. Rev. Lett. 88, 176804 (2002). 
'' J. Lehmann, S. Kohler, P. Hanggi, and A. Nitzan, Phys. Rev. Lett. 88 , 228305 (2002). 

* J. Lehmann, S. Kohler, P. Hanggi, and A. Nitzan, |cond-mat/ 0208404 (2002). 
" J. Lehmann, S. Camalet, S. Kohler, and P. HanggiT ^ysics/ 0205015^ (2002). 

V. Mujica, M. Kemp, and M. A. Ratner, J. Chem. Phys. 101, 6849 (1994). 

S. Datta, Electronic Transport in Mesoscopic Systems (Cambridge University Press, Cambridge, 1995). 

A. Nitzan, J. Phys. Chem. A 105, 2677 (2001). 

E. G. Petrov and P. Hanggi, Phys. Rev. Lett. 86, 2862 (2001). 

J. Lehmann, G.-L. Ingold, and P. Hanggi, Chem. Phys. 281, 199 (2002). 

N. L. Manakov, V. D. Ovsiannikov, and L. P. Rapoport, Phys. Rep. 141, 319 (1986). 

M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 

R. Bliimel, R. Graham, L. Sirko, U. Smilansky, H. Walther, and K. Yamada, Phys. Rev. Lett. 62, 341 (1989). 
S. Kohler, T. Dittrich, and P. Hanggi, Phys. Rev. E 55, 300 (1997). 
J. H. Shirley, Phys. Rev. 138, B979 (1965). 
H. Sambe, Phys. Rev. A 7, 2203 (1973). 

A. G. Fainshtein, N. L. Manakov, and L. P. Rapoport, J. Phys. B 11, 2561 (1978). 
P. Hanggi, in Quantum Transport and Dissipation (Wiley- VCH, Weinheim, 1998). 

R. Bliimel, A. Buchleitner, R. Graham, L. Sirko, U. Smilansky, and H. Waher, Phys. Rev. A 44, 4521 (1991). 

T. Dittrich, B. Oelschlagel, and P. Hanggi, Europhys. Lett. 22, 5 (1993). 

S. Kohler, R. Utermann, P. Hanggi, and T. Dittrich, Phys. Rev. E 58, 7219 (1998). 

P. Hanggi, S. Kohler, and T. Dittrich, in Statistical and Dynamical Aspects of Mesoscopic Systems, Vol. 547 of Lecture Notes 
in Physics, edited by D. Reguera, G. Platero, L. L. Bonilla, and J. M. Rubi' (Springer, Berlin, 2000), pp. 125-157. 
M. A. Ratner, J. Phys. Chem. 94, 4877 (1990). 



